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Detection of Saturnian X-ray emission with XMM-Newton 
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Abstract. The giant planet Saturn was observed by XMM-Newton in September 2002. We present and analyse 
these XMM-Newton observations and compare our findings to the Chandra observations of Saturn. Contamination 
of the XMM-Newton data by optical light is found to be severe in the medium and thin filters, but with the thick 
filter all optical light is sufficiently blocked and the signal observed in the this filter is interpreted as genuine X-ray 
emission, which is found to qualitatively and quantitatively resemble Saturn's Chandra spectrum very well. 
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1. Introduction 

Most of the larger solar system objects are now known 
to emit X-rays via some variety of different X-ray emis- 
sion mechanisms. The most prominent example is the gas 
giant Jupiter, whose X-ray emission is dominated by au- 
roral emission, produced by char ged particles enterin g 
the planet's magnetosphere (e.g.. iMetzeer et all 11983^). 
A recent Chandr a observation of the gas giant Saturn 
l|Ness et al.ll2004 resulted in a definitive detection of X- 
ray emission also from this planet, thus co nfirming a tenta- 
tive ROSAT detection reported earher bv lNess fc Schmitt 
l)200Q|) . however, at a level much lower than observed from 
Jupiter. Further, unlike Jupiter, Saturn's X-ray emission 
is not concentrated in the polar regions, and in fact, the 
detected level of Saturnian X-ray emission is consistent 
with the observed level of Jupiter 's equatorial emission 
ijNess et al.ll2004tlWaite et al.ll996l) . The emission mecha- 
nism consistent with spectral and spatial properties of the 
observed X-ray emission was found to be elastic scattering 
and fluorescent scattering of solar X-rays, however, for this 
to be the c ase, the X-ray alb edo of Saturn has to be unusu- 
ally high. iNess et all l)2004|) estimated the X-ray albedo 
required to explain the measured X-ray flux by scattering 
processes and found a value of > 5.7 ■ 10~^, which is abou t 
a factor 50 higher than for the moon ijSchmitt et alJl99l|) . 
Since models of combinations of scattering processe s for 
Jupiter's equatorial emission jMaurellis et al.ll2000l) un- 
derestimate the observed flux level l|Waite et al.lir99i by 
a factor of 10, the scatter process scenario has to explain 
a high X-ray albedo in both cases, Jupiter's equatorial 
emission and Saturn's total emission. 
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Saturn was also observed by XMM-Newton in 
September 2002. We present and analyse these XMM- 
Newton observations and compare our findings to the 
Chandra observations of Saturn. We concentrate on the 
detection and spectral properties of the X-ray photons, 
while the light curve provides only little information due 
to its short duration. 



2. Observations and Data Analysis 

2.1. Observations 

Saturn was observed with XMM-Newton on September 
10, 2002 for a total of 60ksec. The observations were split 
in three separate parts of almost equal length, with dif- 
ferent fllter settings used for the FN and MOS detectors. 
The MOS detectors were operated with medium filter -|- 
fuU frame, medium filter + large window, and thin filter -1- 
large window. At the position of Saturn no X-ray photons 
can be extracted from the MOS detectors, because pixels 
with high optical load are not read out. The RGS is not 
affected by optical light contamination, but the count rate 
is too low to obtain useful spectra. We therefore use only 
the EFIC-PN detector for our analysis, and the obser- 
vation details are summarized in Table ^ We inspected 
all three FN-observations, but the data taken with the 
thin filter are close to useless and the data taken with the 
medium filter are severely contaminated by optical light. 
However, in the data taken with the thick filter no obvi- 
ous signs of optical contamination are apparent, consistent 
with our expectations about the optical blocking power of 
the thick filter. 

Due to the small apparent motion of Saturn 5", 
compare to the instrument half-power diameter of 15") 
during the 20 ksec observation interval and the high sensi- 
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tivity of the PN-detector we could directly identify Saturn 
on chip jj=A without any need of a transformation fol- 
lowing the apparent planetary motion (see Fig. ^ upper 
panel). After having found emission at Saturn's position 
we carrie d out a tra,nsform ation procedure (described in 
detail by iNess et al]l2004(l . which transforms all recorded 
events into a Saturnocentric coordinate system, and con- 
structed an image in this new coordinate system. We then 
extracted all photons within an extraction radius of 25" 
around the central position of the transformed coordinate 
system, where an enhancement of photons can be imme- 
diately recognized (cf. Fig.^. In a circular detect cell we 
extract 162 photons while from the background (extracted 
from an adjacent 80" x 200" box) we expect only 50.2 pho- 
tons. With a total of 112 ±13 counts we have therefore 
obtained a highly significant detection. 

Table 1. Overvation details for Saturn (only EPIC/PN). 

ObsID 0089370501 

Exp. time 24024 ksec 

Start Time 2002-10-01 10:35 

Stop Time 2002-10-01 17:15 

FN filter THICK FILTER 

on-time (FN) 21047 ksec 

ObsID 0089370601 
Exp. time 
Start Time 
Stop Time 
FN filter 
on-time (FN) 
ObsID 0089370701 
Exp. time 24023 ksec 

Start Time 2002-10-02 00:37 

Stop Time 2002-10-02 07:17 

FN filter THIN FILTERl 

on-time (FN) 20962 ksec 

Angular diam. 
distance (Earth) 
distance (Sun) 
inclination 



24023 ksec 
2002-10-01 17:35 
2002-10-02 00:15 
MEDIUM FILTER 

20966 ksec 



18.1" 
9.2 AU 
9.0 AU 

-26.4° 
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Fig. 1. Photon positions on PN chip (upper panel; 
with expected path of Saturn) and transformed positions 
in "saturnocentric" coordinates (bottom). 162 counts are 
extracted in a circle with radius 25" and 409 counts in an 
80" X 200" background box. 



3. Results 

We analyzed the XMM-Newton EPIC-PN data of Saturn 
in the same fashi on as the Chandra data as described by 
iNess et all l |2004l) . Since the recorded count number is ex- 
actly the same (!), the detection significance in both data 
sets is very similar; note, that the background values dif- 
fer somewhat. Since the angular resolution of the XMM- 
Newton data is lower, we are not able to locate the X-ray 
emission on Saturn's apparent disk (diameter 18.1") from 
our XMM-Newton observations. 

We extracted the light curve with different time bins but 
found no significant variability. The net count rate is 
(5.3 ± 0.6) • 10^'^ cps. Only half a rotation is covered with 
the short observation time and no phase variability can be 
tested. 



Is the signal recorded in the thick filter due to X-rays 
or also due to optical contamination? In order to address 
this issue we analyzed the EPIC-PN medium filter data 
in precisely the same fashion as the thick filter data 
and extracted a background-subtracted spectrum of the 
photons attributed to Saturn (cf. Fig. I^J. The strong 
signal increase towards lower energies is the indicator of 
the severe contamination due to optical light. In contrast, 
carrying out the same procedure with the thick filter data 
results in a spectrum looking totally different (cf. Fig. (S)). 
The thick filter spectrum does not exhibit any increase 
towards lower energies as expected from a genuine X-ray 
spectrum, since the effective areas decrease towards 
lower energies. Further, the thick filter spectrum appears 
very similar to the recorded Chandra spectrum, which is 
overplotted in a light color after scaling to our exposure 
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Fig. 2. Spectrum extracted from the PN observation using 
the medium filter (ObsID 0089370601). Clear indications 
of optical contamination can be identified. 



time and extraction area. The signal is lower due to 
the lower effective areas of Chandra mirrors, with the 
Chandra spectrum appearing somewhat shifted towards 
higher energies. This might be due to some optical 
loading in the Chandra observation, an e ffect that could 
not fully be excluded bv lNess et all pOO^) . A rather weak 
emission feature appears at '^l.SkeV, but is not signifi- 
cant; interestingly it is also seen in the Chandra spectrum. 

Us ing XSPEC we ca rried out spectral modeling simi- 
lar to iNess et all ()2004|) , who found acceptable spectral 
fits with a (physically unmotivated) black body model 
and a combined MEKAL model plus a fluorescent line of 
oxygen. MEKAL contains continuum and line emissivities 
from coUisionally ionized plasma in thermal equilibrium. 
This model is supposed to represent the spectrum of the 
solar corona and the model parameters are the equilib- 
rium temperature and elemental abundances. Given the 
low signal-to-noise of our XMM-Newton data, we can 
only check to what extent the XMM-Newton and Chandra 
spectra are consistent with each other. We rebinned the 
XMM-Newton spectrum to contain at least 15 counts per 
bin, necessary to remain outside the Poissonian regime, 
otherwise a non-standar d statistic al treatment is neces- 
sary re.ef.. rCash.l97fl.: .Ness fc Wichmann 20021 In addi- 
tion we applied the Cash statistics provided by XSPEC 
with the original spectrum and found consistent results. 
In order to present a concise goodness-of-fit parameter we 
here present our results from fits. Our best- fit black 
body model yields a temperature of kT= 0.16 ± 0.03 keV 
(Xred ~ ^-^^ with 9dof), consistent with the temperature 
found from the Chandra observation (O.lSkeV). Instead 
of a MEKAL model we chose an APEC model to describe 
an incident solar spectrum. Assuming solar abundances we 
find a temperature of kT= 0.29 ± 0.05 keV {Xrcd = 0-79 
with 9dof), a little cooler than the temperature found 
from the combined MEK AL/fluorescent line model from 



the Chandra spectrum (kT= 0.39 ± 0.08 keV). A shghtly 
better fit is obtained by introducing an oxygen fluorescent 
line at 527 eV (modeled as a narrow emission line, only 
instrumentally broadened) . With this combined model we 
obtain an APEC temperature of kT= 0.33 ± 0.08 keV, 
consistent with the Chandra results (Xred = '-'•^l with 
8dof). The fit results are summarized in Tableland in 
the last column we list the model fluxes, integrated in 
the wavelength interval 0.1-2 keV. For an overview of the 
available X-ray spectra of Saturn we plot the rebinned 
spectrum with the APEC model (dashed) and the best-fit 
model of the combined APEC/fluorescent line (solid grey) 
in the bottom panel of Fig. |3| 



Table 2. Saturn 


- spectral fits 


and fluxes 






Model 


kT (keV)" 


X?od / dof* 


fiux'= 




EPIC/PN 






(erg/cm.'^ 




black body 


0.16 ± 0.03 


1.09/9 


1.6610" 


-14 


APEC^ 


0.29 ± 0.05 


0.79/9 


1.62-10" 


-14 


APEC^-K 










narrow line'' 


0.33 ± 0.08 


0.41/8 


1.58-^ ■ lO"'" 


Chandra 










black body 


0.18 


0.7/10 


0.44-10" 


-14 


MEKAL+ 


0.39 ± 0.08 


i 


0.55-10" 


-14 


narrow line 




0.9/9 


0.13-10" 


-14 


ROSAT 






1.9-10" 


14 


"90% errors 


''degrees of freedom 






solar abundances 


''at 527 eV, delta profile 


■^0.1- 


2keV 


■'^l.Se-lO-'" (APEC)-F0.22-10-''' 


(fluorescent line) 





While none of the considered spectral models may 
be physically correct, they allow a reasonable accurate 
estimate of the X-ray flux recorded by XMM-Newton. 
An inspection of Table El shows, that an apparent X- 
ray flux between 1.58 — 1.66 ■ 10"^^ ergcm"^ s"^ was 
recorded by XMM-Newton. This compares well with 
the ROSAT flux of 1.9 ■ 10""^"^ ergcm^^ s"^ (in the soft 
ROSAT band 0. 1-0.55 keV) reported bv iNess & Schmitd 
(2000), but is above the level measured by Chandra of 
0.68 - 10^^'' ergcm^^ s^^. Thus an interpretation of the 
signal recorded by XMM-Newton in the thick filter as 
genuine X-ray emission yields flux values consistent with 
earlier observations of X-ray emission from Saturn as well 
as an X-ray spectrum consistent with that recorded by 
Chandra. We therefore conclude that indeed true X-ray 
emission from Saturn has been recorded by XMM-Newton. 



4. Discussion and Conclusions 

We analyzed XMM-Newton observations of Saturn. 
Because of Saturn's visual magnitude of my — 0.9 at the 
time of our observations, the data taken with the thin 
and the medium filters are severely contaminated by opti- 
cal light. However, the data taken with the thick filter are 
exclusively X-ray photons originating from Saturn. Thus, 
X-ray emission from Saturn has been established to be 
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Fig. 3. Extracted spectra for the background and 
source+background (top). The Chandra spectrum, scaled 
to our exposure time and source extraction area, is over- 
plotted. Bottom: Rebinned spectrum containing at least 
15 counts per bin with two models obtained with XSPEC. 
The best- fit model consists of an APEC model and an oxy- 
gen fluorescent line (see Table EJ. 
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significantly weaker than for Jupiter; the reported flux 
levels range from 0.68 • 10^^^ ergcm^^ s^^ from Chandra 
to 1.9 • 10~^'* ergcm~^ s~^ in a marginal ROSAT detection 
l|Ness fc Schmitti2000|) . At least between the Chandra and 
XMM-Newton observations substantial variablity seems 
to have taken place, a fact hardly surprising for almost any 
X-ray source. The spectral models found to be consistent 
with the XMM-Newton data are also consistent (to within 
the errors) with the results from the Chandra observation. 
Any possibly remaining problems with optical loading do 
not seem to significantly affect the results. Further insights 
into Saturn's X-ray production process require substan- 
tially deeper pointings than presently available. 
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